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ABSTRACT 


The theory of anisotropic Heisenberg ferromagnets is developed in the random phase 
and spin wave approximation. It is shown that thermodynamic quantities such as mag- 
netization can be expressed in analytic form for all temperatures. To demonstrate the 
applicability of this technique the magnetization for body -centered, face-centered, and 
simple-cubic lattices was calculated. 
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SUMMARY 

The theory of anisotropic Heisenberg ferromagnets is developed in the random phase 
and spin wave approximation. It is shown that thermodynamic quantities such as mag- 
netization can be expressed in analytic form for all temperatures. To demonstrate the 
applicability of this technique the magnetization for body- centered, face- centered, and 
simple-cubic lattices was calculated. 


INTRODUCTION 

The Ising and Heisenberg models of ferromagnetism have been the subject of inten- 
sive theoretical investigations for many years (ref. 1 and references contained therein). 
Because of the dissimilar behavior of these models, Dalton and Wood (ref. 1) have con- 
sidered an intermediate approach. Using the Green's function method and series expan- 
sion techniques, they studied the critical behavior of their model in both two and three 
dimensions for intermediate values of the anisotropy parameter 77(0 < 77 < 1). 

When treating the thermodynamic quantities of an anisotropic Heisenberg ferromag- 
net, it is usually necessary to use series expansions for the temperature regions of inter- 
est. A purpose of this report is to describe a method of evaluating certain temperature- 
dependent sums or (integrals) that occur in the theory of magnetism without necessitating 
the use of series approximations or extensive computer summations. A single analytical 
expression for the magnetization has been found which is valid for all temperatures. The 
expression, moreover, can easily be extended to include any magnitude of spin. It is 
hoped that the methods applied herein will be found useful in the evaluation of the thermo- 
dynamic quantities in the theory of magnetism and other cooperative phenomena. 
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The method is applied to compute the temperature dependence of the magnetization 
for various values of 77. The simple-cubic (sc), body-centered-cubic (bcc), and face- 
centered- cubic (fee) lattices are considered. 

The anisotropic model of Dalton and Wood is introduced in the random phase approx- 
imation, and an expression for the magnetization is formulated. The calculations for 
the bcc, fee, and sc lattice cases are given. The critical temperature is determined, 
and the results are discussed. A detailed example of the evaluation of the integrals that 
occur in the lattice calculations is given in appendix A. 


ANISOTROPIC FERROMAGNET 

The nearest-neighbor exchange interaction model for the anisotropic ferromagnet in 
the absence of an external field, as described by Dalton and Wood, is specified by the 
Hamiltonian 


H = - 




( 1 ) 


where the sums extend over all sites in the lattice, S j is the spin operator for S = 1/2 
for a spin at site i, and is the exchange integral. (All symbols are defined in appen- 
dix B. ) The Ising and Heisenberg models are obtained by setting 77 = 0 and 1 , respectively. 

Starting with the Green’s function Dalton and Wood derived the correla- 

tion function in the random phase approximation. Their results are 
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The sum in equation (2) goes over all N lattice vectors in the first Brill ouin zone. 
To obtain the spontaneous magnetization we make use of the relation 


s ?g - s(s + » - s g - ( s g) 


For spin 1/2, 




and, hence, the magnetization is given by 


(S z ) +1(1 + 2$) _1 


(4) 


where 



To calculate the magnetization as a function of temperature the sums in equation (5) 
must be evaluated over all values of k in the first Brillouin zone of the appropriate lat- 
tice. Except at the low and high temperature limits, numerical methods are usually 
used. However, it is shown below that these sums can be evaluated analytically over the 
entire temperature range. 

Consider only crystals with cubic symmetry. For these lattices, one usually re- 
places the sum in equation (5) by an integral (ref. 2) 
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where v is the volume per site: sc, v = a , bcc, v = (l/2)a ; fee, v = (l/4)a . 
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For the sc, bcc, and fee lattices, equation (6) can be written in the form 


where for sc, 


for bcc, 


and for fee, 
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x = k x a, y = kya, z = k z a 


(8a) 


x = - k a, y = - k a, z = - k a 
2 x 2 y 2 z 


(8b) 


x =' - k a, y = - k a, z = - k a 
2 x 2 y 2 
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Using the identity 


coth(7rt) = (wt)" 


Equation (5) is written in the form 
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Here r = Rjtt/(S z ) , the reduced temperature t = kT/zJ(0), z is the number of nearest 
neighbors, and c.c. indicates the complex conjugate. 


BODY-CENTERED-CUBIC LATTICE 

Consider first, the case of a bcc lattice, where J(k ) can be written as 


J(0) V x 2 J \ 2 / 


— r z 
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From equation (8b), it follows that equations (10) and (11) have the form 
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With the use of Laplace transforms, the summation in equation (14) can be reduced 
to integrals, and the method is outlined in appendix A. Integrating equations (13) and 
(14) yields the following 
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where K(k) is a complete elliptic integral of the first kind, 


( 16 ) 


(17) 

and Q = <S z )/r. By substituting equations (15) and (16) into equation (9), the magnetiza- 
tion is obtained from equation (4). 

No difficulty arises in the limiting case where 77 = 0, nor does it have to be con- 
sidered separately as is the case when low and high temperature expansions (ref. 1 ) are 
used. By putting 77 = 0, the familiar Ising result for magnetization is obtained. 



<S Z > = i tanh(Q) 
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(18) 


FACE-CENTERED-CUBIC LATTICE 

The treatment for the fee case is analogous to the bcc calculation. The results 

are 
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where 
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SIMPLE-CUBIC LATTICE 


The results for the simple- cubic lattice are obtained in a similar fashion: 
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p _ 3 - 77 cos(x) 
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= cos(k a) + cos(k a) + cos(k a) 
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Figure 1 shows the magnetization <S Z ) for each of the cubic lattices plotted as a 
function of the reduced temperature r. Numerical values were obtained from equa- 
tion (4) for several values of the anisotropy parameter rj. 


CRITICAL REGION 


The temperature at which (S z > = 0 is called the Curie temperature, and its value 
can be obtained from the preceding equations. For example, consider the bcc lattice. 
By substituting equations (15) and (16) into equation (4), one finds 


2r =< 


• fe) K(k) 


-1 


2 3 2 

+ Q coth(Q) - 1 + -2-^- csch 2 (Q)coth(Q) - - 2 - 
8 8 


(27) 


z 3 2 

As (S ) approaches zero, terms such as Q coth(Q) and Q csch (Q)coth(Q) became in- 
determinate. Use of L 'Hospital's rule yields 


q^Q coth(Q) = 1 


(28) 


q^qQ 3 c °th(Q)csch 2 (Q) = 1 


(29) 


and 


2T c ,fe»T (3C 

Equation (30) is the same as obtained by Dalton and Wood by evaluating the sum over 
k-space at the critical temperature in the manner of Watson (ref. 3). Thus, the Curie 
temperature can be determined without the use of any expansions. The same procedure 
can be applied to the other cubic lattices. 

In table I we give the values of the Curie temperature t q for various values of the 
parameter 77 . As expected, we find that our results for 77 = 1 agree with the Curie 
temperature of a Heisenberg ferromagnet calculated (ref. 2) in the random phase ap- 
proximation. 

The preceding calculations were limited to positive values of the anisotropy param- 
eter 77 . If 77 were allowed to become negative, the model would then represent a 
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TABLE I. - REDUCED CURIE TEMPERATURE FOR 


CUBIC LATTICES, NEAREST-NEIGHBOR 
INTERACTIONS FOR SPIN 1/2 


Lattice 

Anisotropic 

Reduced Curie 


parameter, 

temperature, 


V 

kT c /J(0) 

Body- centered- cubic (bcc) 

0 

0. 5000 


.5 

.4830 


1.0 

.3589 

Face- centered- cubic (fee) 

0 

0. 5000 


.5 

.4869 


1.0 

' .3718 

Simple- cubic (sc) 

! 0 

1 

0. 5000 


1 -5 

; . 4776 


1.0 

.3297 




- . 


mixed ferromagnetic and antiferromagnetic system. However, the thermodynamic quan- 
tities are independent of the sign of 77 , in contrast to the anisotropic antiferromagnetic 
two- sublattice model which has been investigated by the present authors (unpublished data 
obtained by L. Flax and J.C. Raich of Lewis). The absence of the mixed state follows 
from the random phase approximation used for decoupling the Green's function equations 
of motion. Other decoupling schemes, such as Callen's (ref. 4), would give different re- 
sults . 


SPIN WAVE APPROXIMATION 

The usual assumption in spin wave theory is that the spin waves are independent of 
each other. This supposition can be expected to be valid only as long as the number of 
reversed spins is small. According to Keffer and Loudon (ref. 3), Dyson (ref. 5), and 
Oguchi (ref. 6 ), the error in the magnetization when the spin wave interactions are ne- 
glected is less than 5 percent at the critical point. 

The nearest- neighbor interaction model of the anisotropic ferromagnet is described 
by equation (1). The spin operators are expressed in the forms obtained by Mallev 
(ref. 7): 
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s t --<3s> 1/2 (v^) 

S+ = (2S) 1//2 b+ 

Sj = -S + b+b f 


where the operators b^ and b^ are to be regarded as the creation and annihilation op- 
erators which satisfy Bose commutation rules. 

The Fourier transforms of bjT and bj are introduced where 
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(32) 


where N is the total number of lattice sites and the sums are over the first Brillouin 
zone. In terms of these operators the transformed Hamiltonian becomes 


H = -NJ(0)S z + 2J(0)S 
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We now construct the Fourier-type Green’s function where there are no interactions. 
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Following reference 4, we need the equation of motion for the single spin wave Green’s 
function ^a^ ; ajj^, where the double brackets « » indicate Fourier transforms of 
the Green’s function, and the single brackets < ) indicate averages over the canonical 
ensemble at temperature T. Hence, the equation of motion for the single spin wave is 


E 



Neglecting the interacting terms, we obtain 


(35) 



where 



= 2SJ(0) 



J(k)“ 

V 

J(0)_ 


(36) 


(37) 


is the energy of a noninteracting spin wave with wave vector k. 

In order to determine the thermodynamic properties, we must solve for the correla- 
tion function, which is given by 


<B;A> = lim 


yO oo 

,lo / (« A ; B >>E«7« A ; B »E-i e ):^r7 

</ oo 6-1 
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for the case t = t*. The result is 



(39) 


The average value of <S Z ) can be obtained from equation (31). 
identical to the Bloch expression. For spin 1/2 it can be written as 


As is expected, it is 


<S Z > 


1 

2 


N 



(40) 
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The solutions for 4* can be obtained from substituting equation (37) into equation (39) 
and using the same procedures as for the random phase approximation. Figure 2 shows 
plots of <S Z > as a function of the reduced temperature for each of the cubic lattices. It 
can be shown that for t]= 1, the Bloch law for magnetization (ref. 8) is obtained 

at low temperatures for the behavior of <S Z ) with temperature. In the extreme aniso- 
tropic case r}= 0, one can obtain the familiar e _1 ^ T law at low temperatures. 

CONCLUDING REMARKS 

An analytical procedure has been outlined for the calculation of the magnetization for 
an anisotropic ferromagnet in the random phase and spin wave approximation over the en- 
tire temperature range of interest. It is hoped that this technique is found useful for 
other applications also. 

The random phase approximation predicts that the spontaneous magnetization shows 
a second-order phase transition independent of the anisotropy. The effect of the aniso- 
tropy parameter is to shift the critical temperature. The critical temperature de- 
creases smoothly as rj increases. Physically this is due to the energy it takes to re- 
verse a spin. For the Ising model, since there is only one component present, an energy 
gap of magnitude 2J(0) is needed to create an excitation. In the isotropic model, the gap 
vanishes due to spin waves. These physical (intutive) arguments are supported by nu- 
merical calculations presented in this report. 

The spin wave approximation gives information regarding the temperature behavior 
of the magnetization for the Heisenberg and Ising models. The results are in agreement 
with the random phase approximation for the appropriate limiting cases. 

A comparison of the analytical technique with the direct numerical summation over 
the first Brillouin zone was made. The results agreed to four places. In addition, prob- 
lems of convergence are avoided by the analytical method. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 6, 1969, 

129-02-05-14-22. 
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APPENDIX A 


EVALUATION OF l 2 

The calculation of the integral Ig involves the summation of an infinite series. 
These summation procedures are discussed by several authors (refs. 9 and 10). Here 
we shall employ some of these techniques and extend them to evaluate Ig. 

As an example of the technique used, consider the face- centered- cubic (fee) lattice. 
The calculations for the other two cubic lattices follow in the same manner. For the fee 
lattice we can write 


I 2 (fcc) = 


3t 


* 2 <s z > 



r Pv K(kj) 

[ Jo (r+?j) + 3ir 


dx + c. c. 


(Al) 


R=1 


k? = 477 


(3 + t\ eos^x) + 3ir 


(3 + v + sir)' 


(A2) 


The complete elliptic integral of the first kind can be represented in series form as 


K(k) = 



V 


< 1 


where 


a = [1 _• 3 • 5. . . (2n - 1)] 
n [2-4*6... (2n)] 


(A3) 


(A4) 


a o = 1 

Only the first two terms of the series (eq. (A3)) are necessary for four-place accuracy. 
There is no difficulty in extending it to as many terms as needed. 

Taking the first two terms of equation (A3) into equation (Al), integrating the com- 
plex function with respect to x, and rearranging terms we find that 
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The following Laplace transforms are used: 


(y 2 + R 2 / 


= ^ t sin (y t ) + yt cos(yt)]| 
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(A12) 



The other sums in Ig and Ig are obtained in a similar way. Integrating equation (A13) 
and its counterparts for the other sums leads to equation (20). 
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APPENDIX B 


SYMBOLS 


bcc body centered cubic 

bp creation operators 

b^, a^ annihilation operator 

excitation energy 
fee face centered cubic 

g, p defined in the text 

H Hamiltonian 

i,j lattice sites 

J(k) defined in the text 

J- exchange interaction 

K(k) elliptic integral of the first kind 

k Boltzmann's constant 

if Laplace transform 

Q <s z >/t 


R sum index 

Sj spin operator at site i 

<S Z > magnetization 

sc simple cubic 

T temperature 

v volume per cit 

y, c defined in the text 

z number of nearest neighbors 

/3 1/kT 

r defined in the text 

r] anisotropy parameter 

r reduced temperature 

(( » Fourier transform of the 

Green's function 

< ) conical average of the operator 
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